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^ ■ Abstract 

C I Following the work of Anily et al., we consider a variant of bin packing, called BIN PACKING 

I ^ ' WITH GENERAL COST STRUCTURES (GCBP) and design an asymptotic fully polynomial time 

approximation scheme (AFPTAS) for this problem. In the classic bin packing problem, a set of 
one-dimensional items is to be assigned to subsets of total size at most 1, that is, to be packed into 
unit sized bins. However, in GCBP, the cost of a bin is not 1 as in classic bin packing, but it is 
a non-decreasing and concave function of the number of items packed in it, where the cost of an 
empty bin is zero. The construction of the AFPTAS requires novel techniques for dealing with 
small items, which are developed in this work. In addition, we develop a fast approximation algo- 
(/3 [ rithm which acts identically for all non-decreasing and concave functions, and has an asymptotic 

^ , ' approximation ratio of 1 .5 for all functions simultaneously. 

1 Introduction 

^ I Classic bin packing 1231 |9l |7l HI is a well studied problem which has numerous applications. In the 

lO ' basic variant of this problem, we are given n items of size in (0, 1] which need to be assigned to unit 

\^ ■ size bins. Each bin may contain items of total size at most 1, and the goal is to minimize the number 

O I of bins used. 

Consider the following possible application. A multiprocessor system, where each bin represents 
one processor, is available for one unit of time. However, a processor that executes a large number 
of short tasks causes the system a larger load than a processor that executes a smaller number of 
^ I long tasks, even if the total duration of the tasks is equal in both cases. This is one motivation to the 

■ problem BIN PACKING PROBLEM WITH GENERAL COST STRUCTURES (GCBP) that we Study here. 

The problem has additional applications in reliability, quality control and cryptography fT]. 

In the problem GCBP, the cost of a bin is not a unit cost, but depends on the number of items 
actually packed into this bin. More precisely, we define the problem as follows. The input con- 
sists of n items / = {1,2, ... ,n} with sizes I > si > S2 > ■ ■ ■ > Sn > 0, and a function 
/ : {0, 1, 2, . . . , n} ^ M.Q, where / is a monotonically non-decreasing concave function, for which 
/(O) = 0. The goal is to partition / into some number of sets Si, ... , Sm, called bins, such that 
X]je5i ^ 1 for 1 < i < m, and so that /(l^'il) is minimized. We say that a function / 
is valid if it has the properties above, and an instance of GCBP is defined not only by its input item 
sizes but also using the function /. We assume that /(I) = 1 (otherwise we can apply scaling to the 
cost function /). 

Anily, Bramel and Simchi-Levi lU introduced GCBP and described the applications in detail. 
We describe their results in what follows. Further results on GCBP appear in lH, but these additional 



o 



'Department of Mathematics, University of Haifa, 31905 Haifa, IsraeL leaSmath . half a . ac . il. 
^Chaya fellow. Faculty of Industrial Engineering and Management, The Technion, Haifa, Israel. 
levinasSie . technion .ac.il. 



1 



results are not related to this paper. A related model was studied by Li and Chen |[T9l . In this model 
the cost of a bin is a concave and monotonically non-decreasing function of the total size of items in 
it. 

For an algorithm A, we denote its cost by A as well. The cost of an optimal algorithm is de- 
noted by OPT. We define the asymptotic approximation ratio of an algorithm A as the infimum 
TZ> 1 such that there exists a constant c, which is independent of the input, so that any input satisfies 
A < TZ ■ OPT + c. The absolute approximation ratio of an algorithm A is the infimum TZ > I such 
that for any input, A < TZ ■ OPT. An asymptotic polynomial time approximation scheme is a family 
of approximation algorithms such that for every e > the family contains a polynomial time algo- 
rithm with an asymptotic approximation ratio of 1 + e. We abbreviate asymptotic polynomial time 
approximation scheme by APTAS (also called an asymptotic PTAS). An asymptotic fully polynomial 
time approximation scheme (AFPTAS) is an APTAS whose time complexity is polynomial not only 
in the input size but also in i. Polynomial time approximation schemes and fully polynomial time 
approximation schemes, which are abbreviated as PTAS and FPTAS, are defined similarly, but are 
required to give an approximation ratio of 1 + e, according to the absolute approximation ratio. 

Anily, Bramel and Simchi-Levi [IJ analyzed the worst case performance of some natural bin- 
packing heuritics when they are applied for GCBP. They showed that many common heuristics for 
bin packing, such as First Fit (FF), Best Fit (BF) and Next Fit (NF), do not have a finite asymptotic 
approximation ratio. Even an application of the first two heuristics on lists of items that are sorted by 
size in a non-increasing order, i.e., the algorithms First Fit Decreasing (FFD) and Best Fit Decreasing 
(bfd), leads to similar results. However, Next Fit Decreasing (NFD) behaves differently, and was 
shown to have an asymptotic approximation ratio of exactly 2. Sorting the items in the opposite 
order gives a better asymptotic approximation ratio of approximately 1.691 (in this case, the three 
algorithms First Fit Increasing (FFl), Best Fit Increasing (BFl) and Next Fit Increasing (NFl) are the 
same algorithm). Note that these heuristics are independent of the specific function /. It is stated 
in fT] that any heuristic that is independent of / has an asymptotic approximation ratio of at least |. 
Therefore, finding an algorithm with a smaller asymptotic approximation ratio, and specifically, an 
asymptotic approximation scheme, requires a strong usage of the specific function /. 

In this paper, we develop an AFPTAS for GCBP. We develop a framework, where the action of 
the scheme for a given non-decreasing concave function / with /(O) = is based on its exact defi- 
nition. We also develop a new approximation algorithm MatchHalf (mh), which acts obliviously 
of /, similarly to the behavior of the algorithms of [IJ. We prove that our algorithm has an asymp- 
totic approximation ratio of at most 1.5 for any non-decreasing concave function / with /(O) = 0, 
improving over the tight bound of approximately 1.691, proved by Anily et al. [1], on the asymptotic 
approximation ratio of NFI. 

The classic bin packing problem is clearly a special case of GCBP as one can set /(O) =0 and 
f{i) = 1 for all i > 1, where the resulting function is monotonically non-decreasing and concave. 
Therefore, GCBP inherits the hardness proof of the classic bin packing problem. That is, GCBP 
cannot be approximated within an absolute factor better than | (unless P = NP). This motivates our 
use of asymptotic approximation ratio as the main analytic tool to study approximation algorithms for 
GCBP. In this metric we design the best possible result (assuming P ^ NP), i.e., an AFPTAS. 

A study of this nature, where approximation schemes are developed for bin packing type problems, 
and in particular, where the complexity of such a problem is completely resolved by designing an 
AFPTAS, is an established direction of research. Studies of similar flavor were widely conducted for 
other variants of bin packing, see e.g. |[T6l[T4ll20ll22l[T3ll . 

Fernandez de la Vega and Lueker lITOl showed that the classic bin packing problem admits an 
APTAS. This seminal work introduced rounding methods which are suitable for bin packing problems. 
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These methods, which were novel at that time, are widely used nowadays. Karmarkar and Karp |[T6l 
employed these methods together with column generation and designed an AFPTAS [16|. In [IT], 
the complexity of two variants of bin packing with unit sized bins are resolved, that is, an AFPTAS 
is designed for each one of them. The first one is Bin packing with cardinality constraints 1 17, 5 1, in 
which an additional constraint on the contents of a bin is introduced. Specifically, there is a parameter 
k which is an upper bound on the number of items that can be packed in one bin. The goal is as 
in classic bin packing, to minimize the number of bins used. The second one is Bin packing with 
rejection |[T2l |3l [TTl, in which each item has a rejection penalty associated with it (in addition to 
the size). Each item has to be either packed or rejected, and the goal is to minimize the sum of the 
following two factors: the number of bins used for the packed items and the total rejection cost of all 
rejected items. Note that prior to the work of 1 13 1, these two problems were already known to admit 
an APTAS |I5l[l2]|3l. The main new tool, used in [13], which allows the design of schemes whose 
running time is polynomial in i, is a treatment for small items using new methods developed in that 
work. The treatment of small enough items for the classic problem is rather simple. Roughly, the 
small items can be put aside while finding a good approximate solution, and can be added later in any 
reasonable fashion. Already in [5], it was shown that if the same treatment is applied to small items 
in the case of cardinality constraints, this leads to poor approximation ratios. Therefore, Caprara, 
Kellerer and Pferschy [5] developed an alternative method for dealing with small items. This method 
still separates the packing of large items from the packing of small items. The scheme enumerates 
a large number of potential packings of the large items, and for each packing, tests the quality of a 
solution that is constructed by adding the small items to the packing in a close to optimal way. The 
enumeration prevents this method from being used for designing algorithms with running time which 
is polynomial in i. The way to overcome this difficulty, used in [13], is to find a good packing of 
large items, that takes into account the existence of small items, and allocates space for them. The 
packing of large items is typically determined by a linear program, therefore, the linear program needs 
to define at least some properties for the packing of small items. Specifically, the linear program does 
not decide on the exact packing of small items, but only on the type of a bin that they should join, 
where a type of a bin is defined according to the size of large items in the bin for bin packing with 
rejection, and on both the size and number of large items, for bin packing with cardinality constraints. 

The problem studied in this paper, GCBP, is more complex than the ones of lfT3l in the sense 
that the cost of a bin is not just 1 . Therefore, even though cardinality constraints are not present, the 
number of items packed into each bin must be controlled, in order to be able to keep track of the cost 
of this bin. In classic bin packing, and other well known variants, forcing all the bins of a solution 
to be completely occupied, results in a perfect solution. To demonstrate the difficulty of GCBP, we 
show the existence of a non-decreasing concave function / with /(O) = 0, for which such a solution 
may still lead to a poor performance with respect to /. 

In our scheme, cardinality constraints are implied by an advanced decision on the cost that needs 
to be paid for a given bin, that becomes a part of the type of the bin. The specific packing of small 
items, which is based on the output of the linear program, needs to be done carefully, so that the 
solution remains feasible, and to avoid large increases in the cost of the solution. An additional new 
ingredient used in our AFPTAS is a pre-processing step, which is performed on small items, where 
some of them are packed in separate bins which are not used for any other items. In typical packing 
problems, bins which contain only very small items are relatively full, and thus the additional cost 
from such bins is close to the total size of these items. However, in our case, such a bin usually 
contains many items, and may result in a high cost. Therefore, our scheme always packs some portion 
of the smallest items separately, before any methods of packing items through a linear program are 
invoked. We show that the increase in the cost of the solution, due to the pre-processing step, is small 
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enough, yet this allows more flexibility in the treatment of other small items, i.e., an additional bin 
would have a small cost compared to OPT. 

The structure of the paper is as follows. In Section |2] we supply examples showing the unique 
nature of the problem GCBP, accompanied with new properties and some properties used in previous 
work. We use all these properties later in the paper. We introduce our fast approximation algorithm 
and analyze it in Section |3] Our main result is given in Section ID 

2 Preliminaries 

In this section we demonstrate the differences between classic bin packing problems, and GCBP. We 
also state some properties proved in lU and HJ to be used later. 

As mentioned in the introduction, common heuristics do not have a finite approximation ratio for 
GCBP LU, and other heuristics have a higher approximation ratio than one would expect. Another 
difference is that sorting items in a non-decreasing order of their sizes is better than a non-increasing 
order. 

A class of (concave and monotonically non-decreasing) functions {fq}qeN that was considered in 
ifn is the following. These are functions that grow linearly (with a slope of 1) up to an integer point 
q, and are constant starting from that point. Specifically, fq (t) = tfort < q and fq (t) = q for t > q. 
It was shown in Ul that focusing on such functions is sufficient when computing upper bounds on 
algorithms that act independently of the function. 

For an integer K > 2, consider inputs consisting of items of two sizes; a = 1 — ^, and b = 

Assume first that there is a single item of size a, and 2K items of size b. NFD packs the large item 
together with K of the small items in one bin, and additional K items in another bin. Consider the 
function fx- The cost of the solution is fxiK + 1) + fxiK) = 2K. A solution that packs all small 
items in one bin and the large item in another bin has a cost of + fK{2K) = K + 1. Thus, even 
though both packings use the same number of bins, the cost of the first packing, which is produced by 
NFD, is larger by a factor that can be made arbitrarily close to 2, than the cost of the second packing. 
Moreover, even though only two bins are used, this proves an asymptotic lower bound of 2 on the 
approximation ratio of NFD (this bound is tight due to LIJ). 

Assume now that there are K items of size a and items of size b. An optimal packing for the 
classic bin packing problem clearly consists of K bins, such that each one is packed with one large 
item and K small items. Using the function fx, this gives a cost of K'^. A different packing collects 
all small items in one bin, and has the cost K ■ fxil) + fxiK'^) = '^K. Since K can be chosen to 
be arbitrarily large, we get that the first packing, which is the unique optimal packing in terms of the 
classic bin packing problem, does not have a finite approximation ratio. Note that this first packing 
would be created by FFD and BFD, and also by FF, BF and NF, if the input is sorted appropriately. 

Throughout the paper, if a specific cost function / is considered, we use OPT to denote the cost of 
an optimal solution OPT for the original input, which is denoted by /, with respect to /. For an input 
J we use OPT (J) to denote both an optimal solution (with respect to /) for the input J (where J is 
typically an adapted input), and its cost. Thus OPT = 0PT(/). For a solution of an algorithm A, we 
denote by m{A) the number of bins in this solution. For an input / we let min(/) to be cost of an 
optimal solution with respect to the function fk for A: = 1, that is, with respect to classic bin packing. 
We let ff^{A{T)) be the cost of an algorithm A on /, calculated with respect to function /fc, and use 
fk{A), if / is clear from the context. 

We further state some lemmas proved in [1] that allow us to simplify our analysis in the next 
section. 
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Lemma 1 [Property 3 in [T]] /i(Nn(/)) = /i(nfd(/)), and therefore w{si) > /i(nfi(/)) - 3. 

iei 

Lemma 2 [Theorem 1 in SJ]J[ Consider a packing heuristic A that does not use information on the 
function f. If the asymptotic approximation ratio of A is at most IZ, for any function fk (for k > 1), 
then the asymptotic approximation ratio of A is at most TZfor any non-decreasing concave function 
f with /(O) = 0. 

A useful packing concept, defined in f\\, is consecutive bins. Recall that we assume si > S2 > 
• • • > Sn- Let Bi, B2, . . . , i?m be the subsets of items packed into the bins created in some solution B 
that packs the items in m bins, where Bi is the i-th bin. The packing has consecutive bins if the union 
Uj<sBj is a suffix of the sequence 1, 2, . . . , n for any 1 < s < m. That is, if the first s bins contain 
n' items, then these are the n' items n — n' + l,...,n — l,n (and thus the smallest n' items). The 
following lemma states that nfi is the "best" heuristic among such with consecutive bins. Consider a 
given input /, the cost function fk and a feasible packing with consecutive bins B. 

Lemma 3 [Corollary 3 in [1][ /fc(NFl(/)) < fk{B{I)). 

A partition of the items (which is not necessarily a valid packing) with consecutive bins is called 
an overflowed packing if for all 1 < i < m, ^ Sj > 1. Clearly, if m > 2, such a packing must be 

infeasible. The following lemma implies a lower bound on the cost of an optimal solution. Consider 
a given input /, a cost function fk, an overflowed packing with consecutive bins B, and a feasible 
packing A. 

Lemma 4 [Corollary 1 in /T/7 fk{B{I)) < fk{A{I)). 

Using these properties, in order to analyze NFI, it is enough to consider the functions fk for k > 1. 
It was shown in [ 1] that the asymptotic approximation ratio of NFI for the function fk (k > 2) is at most 

1 + P The asymptotic approximation ratio of NFI for /i, that is, for classic bin packing, follows from 

00 

the results of [2J and from Lemma[T] This ratio is ^ ^ 1.691. Thus the upper bound of 1.691 

1=1 

ifn follows. In the next section we use these properties to develop a new algorithm. The algorithm 
needs to carefully keep the approximation ratio for /c = 2 while improving the approximation ratio 
for k = 1. 

3 A fast approximation algorithm MH 

In this section we describe a simple and fast algorithm MH, that does not need to know the function / 
in advance. This algorithm is a modification of NFI that tries to combine a part of the relatively large 
items (of size larger than ^) in bins together with one additional item. Note that except for possibly 
one item, NFI packs all such items in dedicated bins. 

As mentioned above, NFI has an asymptotic approximation ratio of at most for the function 
fk with k > 2. Therefore, the difficult case is actually the classic problem. On the other hand, using 
heuristics that perform well for the classic problem, such as FFD, may lead to worse results for > 2 
(which in fact is the case for FFD). Therefore, we define an algorithm that acts identically to NFI, 
except for the usage of a pre-processing step. 
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Algorithm MatchHalf (mh) 

1. Let t be the number of items in / with size in (i, 1] (which are called large items). 

2. Let Mo = { \^] ,---,t}, that is, Mo is the set of smallest [|] large items, and let Mi = 
{!,..., } be the remaining large items. Let S = {t + I, . . . ,n} he called the set of 
small items. 

3. Define the following bipartite graph. One set of vertices consists of the large items of Mq. 
The other set of vertices consists of all small items. An edge (a, b) between vertices of 
items of sizes Sa > ^ and < ^ exists if Sa + < 1, i.e., if these two items can be 
placed in a bin together. If this edge occurs, its cost is defined as c(a, 6) = w{b) (using the 
function w of Section O. 

4. Find a maximum cost matching in the bipartite graph. This matching can actually be found 
using the following greedy process. Insert the items of S into a queue in a sorted order, 
with item t + 1 at the top, and the items Mo are inserted into a queue in a sorted order with 
item t at the top. At each time, let j he the item at the top of the first queue, and i the item 
at the top of the second queue. If Si + Sj < 1, these items are matched, and removed from 
the queues. Otherwise, item j cannot be matched to any item of the second queue (since Sj 
is minimal in that queue), so j is removed from the first queue. This process is done until 
one of the queues is empty, and is performed in Unear time. 

5. Each pair of matched items is removed from /. Every matched pair is packed into a bin 
together. 

6. Pack the remaining items using NFI. 

The greedy process of step|4]finds an optimal matching by a simple exchange argument. We note 
that only (approximately) half of the large items are possibly matched in the pre-processing step. A 
larger fraction may cause an asymptotic approximation ratio above 1.5, as can be seen in the following 
example. Let K he an integer such that K > 2. The input set I consists of K items of size and 
K items of size 1 — Running NFI on this input results in one bin containing K items of size 
and K bins containing one larger item. However, if we match an a fraction (for some < a < 1) 
of the larger items in a pre-processing step, there would be approximately aK bins with two items. 
Consider the function /2. We get /2(nfi(/)) = K + 2, whereas the cost with pre-processing is at 
least aK + K. This would give an approximation ratio of at least 1 + a. 

For the analysis of MH, we use weighting functions. This type of analysis was widely used for 
classic bin packing, and many variants of bin packing. The basic technique was used as early as in 
1971 by Ullman |[23]| (see also |[T5l (1811211 '). We make use of adaptation of the following function 
w : [0, 1] — > M (that is equal to the function Wi{p) defined in f2)\ for any p > 0). We first define 
the well known sequence vTj, i > 1, which often occurs in bin packing. Let tti = 2, and for i > 1, 
TTj+i = TTiiiTi - 1) + 1. Thus 7r2 = 3, TTs = 7, TT^ = 43, ctc. For p e i], we define w{p) = i, 

if /c = TTj — 1 for some i > 1, and otherwise, w{p) = • p. Finally, we let w{0) = 0. Note 
that is a monotonically non-decreasing function. It was shown in [2] that for a given input /, 

^ w{si) > /i(nfd(/)) — 3. Even though both |2| and [1| assume that no zero sized items exist, 

iei 

clearly, the number of bins used by NFD and NFI does not increase as a result of the existence of such 
items, unless all input items are of size zero, and therefore, this property on the weights still holds 
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even if zero sized items are allowed. 

We start with proving the asymptotic approximation ratio for /i . 

Lemma 5 For any input I, m(MH(/)) < | min(I) + 3. 

Proof. We use the following theorem. 

Theorem 6 Consider an algorithm A for classic bin packing. Let wi , W2 be two weight measures 
defined on the input items, Wi : I ^ M, for i = 1,2. Let Wi{I) and W2{I) denote the sum of 
weights of all input items of I, according to wi and W2 respectively, and assume W2{I) < Wi{I). 
Assume that for every input of the algorithm, the number of bins used by the algorithm A is at most 
^2 {!) + T, for a constant value r which is independent of L Denote by Wj the supremum amount 
of weight that can be packed into a bin of the optimal solution, according to measure wi. Then the 
asymptotic approximation ratio of A is no larger than Wj. 

Proof. Given an input / we have A < W2{I) + t. Since an optimal algorithm has OPt(/) bins, with 
a weight of at most Wj in each one of them, we get the upper bound on the weight, according to wi ; 
Wi{I) < Wi • opt(/). Using ^2(1) < Wi{I), we get A < WiOPt{I) + r and the theorem follows. 
■ 

We define a weight measure W2 on items as follows. For every item i, we let W2 (i) = w{si), except 
for small items that are matched to large items in the pre-processing step of MH. These items receive 
a weight of zero according to W2- Let X be the number of bins created by the pre-processing step and 
Y the number of bins created by NFI (i.e., in Step [6] of the algorithm). Let I' be the input after the 
removal of items in the pre-processing step. By Lemma[Tl we have ^ W2{i) = Yl ^(■^i) > — 3. 

On the other hand, every bin created in the pre-processing step has a total weight of 1 , since each such 
bin contains a large item (that has a weight of 1) and a small item of weight 0. Thus ^ W2{i) = X, 

and in total W2{i) > X + Y - 3 = /i(mh) - 3. 

i€l 

Next, we define a weight measure wi. Consider the t large items, and their packing in an optimal 
solution OPT. For any large item a, which is packed in a bin with at least one other (small) item, 
consider the largest small item which is packed with a and denote it by Za- If Za is not well-defined, 
one of the possible items is chosen arbitrarily to be defined as z^. If no such item exists, i.e., a is 
packed as a single item in a bin of OPT, we add an item of size zero to this bin of OPT and define it to 
be Za- Therefore Za exists and is defined uniquely for every large item a. We define the weight of every 
item i as wi{i) = w{si), except for the items Za ior a = 1, ... ,t, for which we let wi{za) = "'^2^°'' . 

In order to show W2{I) < Wi{I), we define a valid matching in the auxiliary graph. This 
matching is based on the packing of OPT. Let Z = {za\l < a < t} and denote a set of the largest [|] 
items in Z = {za\l < a < t}by Z'. We initialize the matching with the items of Z' being matched 
to the large items from their bins in OPT. This matching is valid since by definition of Z, each item 
in this set is packed in OPT in a different bin, with a different large item. If the [|] items matched to 
them are not exactly items , . . . , i, it is possible to replace some large items in the matching by 
smaller large items, until this situation is reached. We have < Si^ for zi G Z \ Z' and Z2 S Z'. 
Since the function w is monotonically non-decreasing, we get ^ w{sza) < 2 Yl ^(■^za)- Let 

ZaeZ Za€Z' 

n 

W{I) = Y wisi). We have W2{I) = W{I) - c{M), where c(M) is the cost of a matching in the 

i=l 

auxihary graph, with a maximum cost, and Wi{I) = W{I) — Y "'^2"°'' > — Y ^i^za) > 

l<a<t ZaGZ' 
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W{I) — c(M) = W2{I), since c(M) is a maximum cost matching on the smallest [|] large items, 
and Yl "^{^za ) is '^ost of one such matching, which we defined above. 

\'-^]<a<t 

Finally, we need to find an upper bound on the total weight in a bin of OPT, according to wi. We 
first consider bins that do not contain a large item. For any item i of size Si = P G (0, we have 
wi{i) < For items of size the weight is 0. Therefore, the total weight of items in such a bin is 
no larger than 1.5 (a tighter upper bound of 1.423 is proved in [2J). 

Consider next a bin which contains a large item. Let a be the large item of this bin, and Za is 
chosen as above. If s^^ = 0, then the only item in the bin that has a non-zero weight according to wi 
is a, and thus the total weight is 1. Otherwise, let j be such that Sz^ G (jTT' j]- '^'^^ other item i in 
the bin (except for a and z^) satisfies W2{i) < '^Sj (since Sj < s^a < j)- If J = ^Tj — 1 for some 
i > 1, we have W2{za) = ^- Otherwise, W2iza) = ^f^za- 

We have a total weight of at most 1 + W2{za) + ^^{l- Sa- s^J < 1 + -^2(^0) + "^(^ - ^z^), 
since Sa > ^. In the first case we use > j^, and get at most 1 + ^ + '^jj- ~ j = i- Iii second 
case we get at most 1 + ^j-Sza, + ^j- — ^'J'^za = ~ ^^^za- Using the same property we get 
at most I again. ■ 

Next, we perform an analysis for functions fk with k > 2. Let / be the original input on which 
MH is executed. Let / denote an input in which every small item, which is matched with a large item 
in the pre-processing step of MH, is replaced with an item of size si. Thus, at most [|] items are 
increased to the size si. We consider the following solutions and compare their costs. The cost of 
the solution of MH on /, with respect to fk, is denoted by Ak{I). The cost of the solution of NFI 
on /, with respect to f^, is denoted by NFlfe(/). The next solution that we consider is an overflowed 
solution that is created for / as follows. The items are sorted by size in a non-decreasing order (that 
is, order by indices in a decreasing order). At each time, a minimum prefix of the items of total size 
larger than 1 is assigned to the next bin. The cost of this solution with respect to fk is denoted by 
Ok{I)- The cost of an optimal solution for /, with respect to fk, is denoted by OPTfc(/). Finally, we 
consider a solution for / with consecutive bins, which is constructed from the overflowed solution for 
I as follows (the construction is similar to the one in HI, except for the treatment of items in /, and the 
fact that the corresponding items in I are simply removed). For every bin of the overflowed solution, 
if the total size of items exceeds 1 (this is the case with all bins except for possibly the last bin, or bins 
with removed items), remove the last item and open a new bin for it. The additional large items of 
j, which existed as smaller items in / and were removed from I, are assigned to dedicated bins. The 
cost of this solution, with respect to fk is denoted by Ck{I)- 

By Lemma [3l we have NFlfc(/) < Cfc(/). By Lemma IH we have Ok{I) < OPTfc(/). We next 
prove two lemmas after which we will be able to conclude Ak{I) < |oPTfc(/) + 3.5. 

Lemma 7 Ak{I) < NFifc(/) + 1. 

Proof. Since all small items of I that are packed in the pre-processing step of MH are large in /, the 
small items packed by NFI in the two algorithms are the same ones, and bins created by NFI in the two 
algorithms are identical, except for bins that contain a large item. If any of the two applications of NFI 
outputs a bin that contains a large item together with other items, we adapt the solution by moving this 
item into a separate bin, this modification cannot decrease the cost of a solution, but it may increase 
the cost by at most 1. The small items bins, resulting from running NFI in both solutions (the solution 
of MH and the solution of NFI, possibly with the modification) are now identical. The remaining items 
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are packed in both solutions either in singles or in pairs. Thus the costs of such bins are equal in both 
solutions (since k > 2). Therefore, the claim follows. ■ 

Lemma 8 Cfe(/) < |Ofc(/) + 2.5. 

Proof. We first modify both solutions so that none of them combines large items with some small 
item in one bin (but the overflowed solution may still have bins with two large items, which are not 
modified here). For the overflowed solution, this may require moving one or two large items from a 
shared bin to a dedicated bin, so it may increase the cost by at most 2. For the other solution, this may 
involve moving one large item to a dedicated bin, and cannot decrease the cost of the solution. We 
consider first the bins with small items, that contain at least A; + 1 items in the overflowed solution. 
For every such bin, its cost is at least k. As a result of moving the last item to a dedicated bin (in the 
process of creation of the feasible solution), an additional cost of at most 1 is incurred. Thus the cost 
increases by at most a factor of |. For any bin containing at most k items, there is no additional cost 
from this step. Note that all bins with large items are in this situation. The cost of bins with large 
items in the overflowed solution with the modification is simply t, no matter how they are exactly 
packed, so packing each one in a dedicated bin does not change the cost. Together with the additional 
[|] large items, the cost of large items becomes [y] < y + ^. Removing small items that do not 
exist in / may only decrease the cost. This proves the claim. ■ 
Using Lemma m we have proved the following. 

Theorem 9 The asymptotic approximation ratio o/MH is at most 1.5. for any non-decreasing concave 
function f with /(O) = 0. 

We have shown above that for = 2 (and similarly, for any constant k), the bound 1.5 is tight. 
Note that the bound 1.5 is tight for /c = 1 as well. Consider an input with N large items of size 
i + 27^, and N{K — 1) small items of size ^ (for large enough N, K, such that N is divisible by 

. r^x AT , ■ ■ . , . ,, ■ N{K-\:)~^- N(K-^) , . ., 

4a). mh creates ^ bms with one large and one small item, — ~ — 2K "^^^ ^^^'^ 

small items each, and ^ bins with one large item. This gives a total cost of + . An optimal 

solution combines K — \ small items with every large item, for a cost of A^. For large enough K, the 
ratio is arbitrarily close to 1.5. It can be seen that this ratio is achieved for any fraction < a < 1 of 
large items that participate in the pre-processing step. 



4 An AFPTAS for GCB P 

In this section we present our main result, that is, an AFPTAS for GCB P. We give a sketch which 
presents the main ideas and technical difficulties, and give the full description of the AFPTAS and its 
analysis later. We first present an auxiliary algorithm called Fractional Next-Fit Increasing. 



4.1 The analysis of FNFI 

We prove a property which is helpful in the design of our AFPTAS. It is related to the property on 
NFI in Lemma [3l but it is stronger since it is proved for any non-decreasing concave function / with 
/(O) = 0, for fractional packing of items. A packing is fractional if items can be cut into pieces, where 
pieces of one item can possibly be packed in different bins. We assume without loss of generality that 
in every fractional packing, every bin contains at most one part of each item. If this property does not 
hold, it is possible to unite parts of items within a bin without changing the cost. 
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We consider an algorithm which creates a fractional packing of the items according to the variant 
of the NFI heuristic, called FRACTIONAL NFI (FNFI). This algorithm sorts items by size in non- 
decreasing order. At each time, a bin is filled completely, before moving on to the next bin. For this, 
we allow the splitting of items into several parts, that is, the last item that is packed in a bin is possibly 
just a part of an item. Consequently, the first item packed in the next bin may be the remaining part 
of the same item. Note that each bin in the output of FNFI contains at most two split items and that in 
total only at most m — 1 items are spht (where m is the number of bins used by FNFI). 

Note that there is no advantage in packing fractions of size zero of items, except for zero sized 
items, which we assume that are split between bins. If a part of size a of an item of size /? > is 
packed in a given bin, we say that the fraction of this item that is packed in this bin is ^. If an item is 
packed in a bin completely, we say that its fraction packed in the bin is 1. The number of items in a 
bin which is packed fractionally is the sum of fractions in it. This number is not necessarily an integer 
and it is unrelated to sizes of these fractional items, but only to their fractions. 

To be able to analyze fractional packings, we next define / for any (real and not necessarily 
integral) value q G [0,n] as follows. We define f{q), for i < g < i + 1, to be (i + 1 — g) • /(i) + 
{q ~ i) ' f{i + !)• The values of / for integer values of q are unchanged. We let f{x) = f{n) for 
any x > n. This function is piecewise linear and continuous, and since it is an extension of a non- 
decreasing concave function on integers, it is monotonically non-decreasing and concave in [0, n]. The 
cost of a fractional packing is calculated according to the generalized function /, using the numbers 
of items packed into the bins as defined above. 

A simple property of FNFI is that it creates bins that are sorted in a non-increasing order of the 
number of items in them. This holds since given two bins ii < 12, bin ii is completely occupied, and 
every item that has a part packed in bin ii has a size no larger than any item that has a part packed in 
bin 12. 

For any non-decreasing concave function / with /(O) = 0, the following lemma states that fnfi 
is the "best" heuristic among packings with fractionally packed bins. Consider a given input /, a cost 
function / and a fractional packing, B. 

Lemma 10 /(fnfi(7)) < f{B{I)). 

Proof. Assume by contradiction that for an input /, a fractional packing B and a function /, we have 
/(fnfi(/)) > f{B{I)). Assume that the bins of B are sorted according to a non-increasing numbers 
of items. If the packing B that satisfies the condition is not unique, consider such a packing B which 
maximizes the suffix of bins that are packed identically to the packing of FNFI. Consider the first bin i 
of B that is packed differently from the packing of fnfi. If bin i is the very last bin of the packing B, 
then the bins 1, . . . , i — 1 are packed as in the packing of FNFI, and therefore, bin i also has the same 
contents for B as it has for fnfi. Therefore we assume that i is not the last bin of B. 

Let J, J ■ + 1, . . . , j' be the indices of items that FNFI packs in bin i (the first and last items, which 
have the indices j' and j respectively, may be packed fractionally in this bin). Let j < ji < j' be an 
index of an item such that B packs a smaller part of ji (possibly of size zero) in bin i than FNFI does. 
Such an item must exist by the following argument. If FNFI fills bin i completely, then since bin i of 
B is packed differently, it cannot have at least the same fraction of every item. Otherwise, FNFI packs 
all the remaining items in bin i, so a different packing of bin i means that some item has a smaller 
fraction in B. 

We next consider the case that there exists an item j2 for which B packs a larger part in bin i than 
the packing of FNFI. Since the two algorithms pack bins 1, . . . ,i — 1 identically, only the items of 
index up to / are available for packing in bins i,i + 1,. . ., where the item of index / may akeady 
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be fractional. Out of these items, FNFI packs a maximum prefix into bin i, so this item must satisfy 
j2 < j- We get that j2 < J < Ji- Since ji ^ j2 by their definitions, we get j2 < ji- 

Denote the fractions of ji and j2 in bin i of by 71 and 72, and the fractions of ji and j2 in bin i 
of FNFI by 61 and 62- We have 5i > 71 > and 72 > (^2 > 0. Since 71 < 61, and bins 1, . . . , i — 1 
are packed identically in both algorithms, there exists a further bin i' that contains a part of item ji in 
the packing of B. Let ei > be the fraction of ji in bin i' of B. 

We would like to swap parts of items in the packing of B, specifically, a part of item ji from bin 
i' with a part of item j2 in bin i. We use jj, to denote the size of the swapped part. There are three 
restrictions on jj.. The resulting fraction of ji in bin i of B cannot exceed the fraction of this item in 
bin i of FNFI, thus /i < (5i — 7i)sji. We can swap at most a fraction ei of ji. Moreover, we can swap 
at most a fraction of 72 — 62 of j2, in order to keep a fraction of j2 in bin i that is at least as large as 
the one in bin i of FNFI. Therefore, we let /x = min{(7i — (5i)sj^, (72 — S2)sj2,£iSj-^}. We adapt B 
by swapping a part of size fi of item ji from bin i' with a part of size fj, from j2 in bin i. By definition 
of all variables, fi > 0, and thus some change occurred. 

Let n-i and n^/ be the original numbers of items in bins i and i' of B. By our assumption > rij/. 
Let Qi and 02 be the fractions of items ji and j2 that are swapped. Since /i = ai • Sj^ = Q2 • sj^, 
and < Sj2, we have ai > 02. Thus, the change in the cost is /(rij — 02 + ai) + f{ni' — ai + 
0:2) — /(^i) + f{ni') < 0, by concavity. As a result of this process, the total number of items in bin i 
remains no smaller than the numbers of items in each of the bins i + l,i + 2, . . .. 

If an item j2 does not exist, it means that bin i has a total size of items that is smaller than the total 
size of items in bin i of FNFI. In particular, it means that bin i is not fully packed. We define 71, 61, 
i' and ei as before. In this case we can define fi = min{(7i — 6l)sj-^^,elSj-^}. We define ai, and 
rij/ as before. Thus, the change in the cost is /(rij + ai) + /(nj/ — ai) — /(nj) + /(nj/) < 0, by 
concavity. 

It is possible to perform this process on bin i multiple times, until there is no item that has an item 
for which a smaller fraction of it is packed in bin i of B than it is packed in the same bin for FNFI. At 
this time these bins become identically packed. 

We next show that this situation, where no item ji exists, is reached after a finite number of swaps. 
For every item ji , it can be performed for every item j2 and for every successive bin. This gives a 
total of at most swaps, and possibly movements of items to bin i without swaps. 

After we reach the situation where bin i is identical for B and FNFI, the bins 1, . . . , z of ;B are 
sorted by a non-increasing number of items. Each remaining bin of B has a number of items that is 
no larger than bin i. Moreover, bins i + 1, i + 2, . . . can be sorted so that the list of bins becomes 
sorted as required. The changes above can only decrease the cost of the solution, and therefore we get 
a contradiction to our assumption. ■ 

4.2 The sketch of the scheme 

We define an item to be a small item if its size is smaller than e and otherwise it is a large item. Denote 
by S the set of small items and by L the set of large items. Our first step is to apply linear grouping 
ifTOl of the large items, that is we sort them by size and we partition them into (almost) equal-sized 
sets of consecutive items (in the sorted list). We pack each item of the set of the largest items in its 
own bin, and we round-up the size of the items in each other set to the largest size of an item in its set. 

We next partition the items in S into S' U S" where S" contains the smallest items such that the 
total size of the items in S" is close to a constant which we define depending on e. The items of 
S" are packed nearly optimally using the FNFI heuristic and packing any split item using a dedicated 
bin. These bins will enable us to use a constant number of bins with an arbitrary content (of items 
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in L U 5') while paying at most e times the cost of the bins which are used to pack the items in S" . 
We note that packing S" using the NFI heuristic is also possible and leads to a similar performance 
guarantee. However, the analysis of using FNFI is simpler. 

Our next step is to approximate the cost function / using a staircase (step) function with 0(log f{n)) 
steps. We use concavity of / to show that this number of steps in the function is sufficient to get a 
1 + e approximation of /. 

We next move on to finding a packing of the items in L U S" (neglecting the largest items which 
are packed in dedicated bins). In such an instance, the linear program, which we construct, allows 
the small items of S' to be packed fractionally. To construct this linear program we define a set of 
configurations of large items (this is the standard definition), and a set of extended configurations 
which also define the space and cardinality of small items in a configuration (this is a non-standard 
idea). The linear program will decide how many bins with a given extended configuration to open 
and what type of bins each small item need to be packed in. These types are called windows, and 
we define them as the pair consisting of the total space for the small items and the total cardinality 
of small items in a bin with this window. Hence in this linear program we have a constraint for each 
size of large items (a constant number of constraints) a constraint of each small item (a linear number 
of such constraints), and two constraints for each type of windows. We apply the column generation 
technique of Karmarkar and Karp [16] to solve approximately the resulting linear program (we use a 
separation oracle which applies an FPTAS for the Knapsack problem with cardinality constraint given 
by L6J). 

Unfortunately the number of fractional entries in a basic solution for this linear program (as we 
can assume our solution is indeed a basic solution), is linear in the number of windows types (plus 
a constant). The number of windows is indeed polynomial in the input size allowing us to solve the 
linear program, but it is not a constant, and we will incur a too large error if we would like to round 
up the fractional solution. 

Hence, we define a restricted set of windows types with a much smaller set of windows, and 
we show how to project cleverly our solution to a new solution which is not worse than the original 
solution, whose support uses only windows from this restricted set of windows. Therefore, when we 
count the number of constraints, we can eliminate the constraints corresponding to windows which 
do not belong to the restricted set of windows. Thus the new bound on the number of fractional 
components in the projected solution is now much smaller. That is, our projected solution which is 
an approximated solution to the original linear program is also an approximated solution to the linear 
program with additional constraints setting the variables to zero if the corresponding window does not 
belong to the restricted set of windows. 

The next step is to round up the resulting projected solution. If a small item is packed fractionally, 
then we pack it in its own dedicated bins. If the fractional solution needs to pack fractional copies of 
bins with a given extended configuration, then we round up the number of such bins. The large items 
clearly can be packed in these bins according to the configurations of the large items. The small items 
are now assigned to windows (by an integral assignment), and not to specific bins. Therefore, our last 
steps are devoted to packing the small items. 

We first place the small items which are packed in a common window type into the bins with 
this window as part of their extended configuration in a round-robin fashion where the small items 
are sorted according to their size (this ensures us that the number of items in each such bin will be 
approximately the same, and the total size of these items in such bins will be approximately the same). 
Hence, the excess of volume of small items in a bin is relatively small (with respect to the total size of 
small items in this bin). In fact it is at most one excess item per bin plus a small volume of additional 
small items (this small volume is due to a rounding we have done when we define the set of windows). 
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The excess items are packed in dedicated bins such that ^ excess items are packed in each dedicated 
bin. The small volume items are packed again in dedicated bins such that these items from | bins are 
packed into one common dedicated bin. The items which are removed from a bin after the process 
of the round-robin allocation are the largest small items of this given excess volume. The resulting 
scheme is an AFPTAS for GCBP, as claimed by the following theorem. 

Theorem 11 The above scheme is an AFPTAS for GCBP. 

4.3 A detailed description and analysis of the AFPTAS for GCBP 

Let < £ < I be such that ^ is an integer. Recall that /(O) = and /(I) = 1. 

The input for this problem includes in addition to the list of items, also the function /. Therefore, 
the running time needs to be polynomial in the following four parameters: n, ^, and the binary 
representations of the numbers in the input, including the item sizes, and the values of / on the 
integers 1, . . . , n. The length of the representation of / is at least log /(n). 

If n < we pack each item into a separate bin. In this case, the cost of the solution is at most 
< (1 + £)OPT + -. We therefore assume that n> -. 

Linear grouping. An item j is large if Sj > e. All other items are small. We denote by L the set 
of large items, and by S the set of small items. We perform linear grouping of the large items. That 
is, if l-Z^I > ^, then for m = ^ we partition Linto m classes Li, . . . ,Ljn such that [|i|e^] = |ii| > 
\L2\ > ■ ■ ■> \Lm\ = [\L\e'^\, and Lp receives the largest items from L \ [Li U • • • U The 
two conditions uniquely define the allocation of items into classes up to the allocation of equal size 
items. For every j = 2, 3, . . . , m we round up the size of the elements of Lj to the largest size of an 
element of Lj. For an item i, we denote by s[ the rounded-up size of the item. If |L| < ^, then each 
large item has its own set Li such that Li is an empty set, and for a large item j we let s'^ = sj (i.e., 
we do not apply rounding in this case). In both cases we have < 2e^|L|. 

For items in Li, we do not round the sizes, and we denote s'j = sj for all j G Li. For j G S 
we also let s'j = sj. We denote hy L' = L \ Li. We consider the instance /' consisting of the items 
in L' U S with the (rounded-up) sizes s'. Then, using the standard arguments of linear grouping we 
conclude opt(/') < opt(/). The items in Li are packed each in a separate bin. We next describe the 
packing of the items in I'. 

Dealing with the set of the smallest items. We define a partition of the set S into two parts S' 
and S", such that S" is a suffix of the list of input items (i.e., a set of smallest items). Specifically, 
if i G 5' and j G S", then s'^ > s'j. Let S" be a maximum suffix {p, . . . , n}, such that S" C S, for 
which the total size is at most 1 + h{s), where h{e) is a function of £ that we will define later. This 
function is defined such that h{£) > ^ is an integer for any valid choice of e. Note that if the total size 
of the small items is smaller than 1 -|- h{e) then we let S" = S and S' = 0. We will pack the items 
from S" independently from other items. That is, there are no mixed bins containing as items from 
S" as items not from S". 

The first packing step of the algorithm is to pack the items of S" using the following heuristic. We 
apply FNFI (processing the items in an order which is reverse to their order in the input). This results 
in 1 + hie) bins, unless 5"' = S. Afterwards, a new dedicated bin is used for every item that was split 
between two bins by fnfi. There are at most /i(e) such items. 

In order to focus on solutions that pack the items of S" as we do, we next bound the cost of a 
solution that packs the items in S" in this exact way (packed by fnfi in separate bins, where split 
items are moved to an additional bin). On the other hand, we relax our requirements of a solution and 
allow fractional packing of the items in S'. The solution clearly needs to pack the items in L' as well 
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(no fractional packing can be allowed for large items). We denote the optimal cost of such a solution 
by OPt'(/'). The motivation for allowing fractional packing of the items of S' is that our goal is to 
bound the cost of solutions to a Unear program that we introduce later, and this linear program allows 
fractional packing of small items that are considered by it, which are exactly the items of S' (while 
the items of S" remain packed as defined above). 

Lemma 12 opt'(I') < (1 + e)OPT(/') + (3/i(e) + 3) • /(i) < (1 + e)OPT + {3h{e) + 3) • /(i). 

Proof. Consider an optimal solution OPT to the following relaxation GCBP'^ of our packing prob- 
lem. We need to pack the items of /' (with rounded-up sizes) but all the items of S can be packed 
fractionally. The difference with the packing OPt'(/') is that items of S" can be packed in an arbi- 
trary way, and not necessarily into dedicated bins, as is described above. In particular, they can be 
packed fractionally. The difference with the packing OPt(/') is the possibiUty to pack the small items 
fractionally. The cost of OPT is clearly at most OPt(/') < OPT. 

We sort the bins of OPT in a non-increasing order, according to the number of items (i.e., the sum 
of fractions of items) packed in the bin (including large items). Let (jj be the total free space in bin i 
that is left after packing its large items in it. This is the space which is used by small items, together 

i 

with all the free space, if exists. Let Sj = ^ (jj. help = min{z|Sj > ^ s'A. The integer p must 

j=i jeS" 

exist since all items of S" must be packed. 

We show that without loss of generaUty, we can assume that all items of S" are packed in bins 
1 , 2 , . . . , p in OPT. To show this, consider an optimal solution to GC B P ' ^ that minimizes the following 
function (among all optimal solutions): the number of existing quadruples (ai, ii, 02, ^2), where ai < 
p < 02, n G <S", i2 G S", and there is a non-zero fraction of item ij packed in bin aj, for j = 1, 2. 
Assume by contradiction that such a quadruple (ai, ii, 02, 12) exists. Let 7 be the fraction of ii in bin 
ai and 6 the fraction of 12 in bin 02. 

Let LL = min{7 • Sj, , 5 • Sj,}. Denote the fractions of ii and 12 of size M by 7' = -7^ and S' = 
We swap a part of size /x of item i2 in bin 02 with a part of size /x of item ii in bin ai. Since Si^ > Si^ 
(recall that S" contains the smallest items), we get that the fractions satisfy 7' < 5' . The number of 
items in bin ai was changed by (5' — 7', and in bin 02 it was changed by 7' — (5'. The sorted order of bins 
may have changed as a result, but bin ai can be moved to an earlier spot while 02 may be moved to a 
later spot, so the set of the first p bins does not change. Moreover, we destroyed at least one quadruple, 
and did not create new ones, since no parts of items of S' were moved to bins 1, . . . , p and no items of 
S" were moved to bins p+l,p+2,.... Let rii and n2 be the numbers of items in bins oi and 02 before 
the change. The change in the cost function is /(rii + — 7') + /(n2 — (5' — 7')) — /(^^i) — /(ra2) > 0, 
since ni > n2, 5' — 7' > 0, and by concavity. Therefore, the resulting solution has a cost of at most 
OPT, and the minimaUty is contradicted. 

If no such quadruple exists then there are two cases. If all bins p + l,p + 2, . . . contain only 
fractions of items of S' (possibly in addition to large items), then all items of S" are in bins 1, . . . ,p 
and our assumption holds. Otherwise, we have that all bins 1,. . . ,p contain no fractions of items in 
S'. In this case, if there are items of S" in any of the bins p + l,p + 2,..., then there must be empty 
space in bins 1, . . . ,p. Parts of items of S" can be repeatedly moved to these bins, until no parts of 
items of S" exist in bins p+l,p + 2,.... In each such step, the number of items in some bin in 
I,. . . ,p increases, and the number of items in some bin mp + l,p + 2, . . . decreases. Sorting the bins 
again after every such step (according to a non-increasing numbers of items) will contain the same 
set of bins in the prefix of p bins, and our assumption holds as well. Due to concavity, and since the 
target bin cannot contain less items than the source bin, every such step cannot increase the cost. 
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We next adapt OPT by creating at most h{e) + 2 additional bins, and move the small items of the 
first p bins into these bins using FNFI (that is, the list of items is processed in a reverse order from 
their order in the input and packed fractionally into bins). Note that this set of small items may contain 
items of S' of total size at most 1 (out of these items of 5", at most one is split between two bins), and 
the total size of items of S" is at most h{e) + 1. We denote this set of items that is moved by S. We 
compute the change in the cost and afterwards adapt the solution further so that it complies with the 
requirement that the items of S" are packed integrally in separate bins, as is done above. 

We define an auxiliary monotonically non-decreasing concave function / as follows. f{x) = 
f{x + ^) — /(^)- Note that /(O) = 0. Consider the p bins of OPT from which the small items are 
removed. Let and Oj denote the numbers of large and small items in these original bins. Clearly, 

< p By removing the small items, the cost of such a bin decreases by f {ai + Vi) — f (ri) > 
f{ai + 7) — /(^) = f{o,i), where the inequality is due to concavity. For every bin which is created 
for small items, if it contains hi small items, its cost is f{bi) < /(6j + = + f{\), where the 
inequality is due to monotonicity. 

Consider now the packing of the items S that is implied by the solution OPT, with respect to the 
function /, and neglecting the large items. The cost of this packing for bin i is f{ai). Let A denote 
the total cost of all the bins that contain items of S, that is, of the first p bins. Let B denote the total 
cost with respect to / of all the bins that are created by FNFI for S. In this case the cost of a bin i is 

P ^ , /i(e)+2 

f{bi). That is, A = ^ /(fli) and B = Yl fih)- By Lemma [TOl (that holds even though the value 

i=l i=l 

/(I) can be arbitrary), we have A> B. 

h{e)+2 p 

Let A denote the difference in the cost for the items of S. We have A = Yl f{bi) — {fi'^i + 

i=l i=l 

ai) - f[n)) < Y ifibi) + /(i)) - E Rai) < [h{e) + 2)f{\) (by the previous claims and 

i=l i=l 

A > B). 

We next convert the packing of small items as follows. If there exists a mixed bin, that is, a bin 
containing items from both S" and S' , we split it into two bins, so that the two subsets of S' and of 
S" are separated. If a mixed bin indeed exists, S' 7^ 0, and the total size of the S" items is more than 
h{e), but not more than h{e) + 1. Therefore, the split bin appears as the h{e) + 1-th bin created by 
FNFI. Moreover, the number of items in the h{e) + 1-th bin is no larger than the number of items in 
every earlier bin. Therefore, if the number of items in the + 1-th bin is N, then the current cost is 
at least f{N){h{e) + 1) and as a result of the split, the cost increases by an additive factor of at most 
f{N). So the multiplicative factor of the increase in the cost is at most 1 + ^^^-^^^-^ < 1 + e where the 

inequality holds by h{e) > ^. 

For a pair of consecutive bins created by FNFI (excluding the bins with items of 5* n S"), if an item 
was split between the two bins, it is removed from these bins and packed completely in a new bin 
dedicated to it. There are at most h{s) such items so this increases the cost by at most h{e) • /(I). At 
this time, the items of S" are packed exactly as in OPt'(/'). 

The total cost is at most (1 + e)(OFf + (/i(e) + 2)/(i)) + h{e) < (1 + e)(OFf ) + ((2 + e)h{e) + 
2 + 2e)/(i) < (1 + e)(OFf) + {3h{e) + 3)/(i) (using e < i). ■ 

We next need to pack the items in /" = I' \ S". Let 6 = min s •. Clearly, for any i G 5" we have 

s'i < 5. Let A = f 

We next consider the instance /". In the temporary solutions, we allow fractional packing of the 
items of S' and we use OPt(/") to denote an optimal packing of /" where small items may be packed 
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fractionally. This does not change the fact that any bin, packed with items of a total size of at most 1, 
can contain a total number of items of at most A even if it contains fractions of items. 

We denote the cost of the bins packed with the items of S" by F{S"). By definition we have 
OPt'(/') = OFt{I") + F{S"). The items of S", if packed by FNFI (which by LemmafTOlis a minimum 
cost packing for them) require at least h{e) full bins, with at least A items in each. Therefore, we 
have F{S") > h{e) ■ /(A). On the other hand, at this time, any other valid bin can contain a total 
number of items of at most A. 

These properties are true unless S' = 0. In that case, only large items remain to be packed, so the 
number of items in any additional bin is at most K In this case we let A = i. 

Approximating the cost function /. Given the function / we compute a staircase function, which 
is an (1 + e)-approximation of /, with 0(logi^j /("^)) breakpoints. That is, we find a sequence 
of integers = A;o<A;i = l<---<A;i=-< ki,^ < ■ ■ ■ < ki = n such that for all 

i = - + 1, . . . ,£ — 1,-we have /(fcj+i) < {l + e)f{ki). The sequence is constructed as follows. We 
define kj = j for j = 0,1, . . . , Every subsequent value kj^i for j > Ms defined as the maximum 
integer t > kj such that f{t) < {l + e)f{kj). Note that this definition is valid since for j > ^ we have 
f{j + 1) < /((I + < (1 + £)f{j), where the first inequality holds by the monotonicity of /, and 
the second inequality holds by the concavity of /. Then, by the definition of the sequence, for every 
i = i i + 1, — 2, we have f(ki+2) > (1 + £)f{ki)- Note that by the definition of this sequence, 
we have i = + logi+e /('^)) ^^'^ £ < n. Let p/\ be such that kp^ > A and kp^-i < A. If 
S' = 0, we have A = -, so kp^ = -. The staircase function, which is an (1 + e) -approximation of 
/, is defined as the value of / for values ki, and it remains constant between these points. 

Constructing tlie linear program. Given the instance /", we let a configuration of a bin C be a 
(possibly empty) set of items of L' whose total (rounded-up) size is at most 1. We denote the set of 
all configurations by C. For each configuration C we define p/\ + 1 < i + 1 extended configurations 
{C, ko), {C, ki), . . . (C, kp^). A bin packed according to an extended configuration (C, kp) has large 
items according to configuration C, and at most kp items in total (that are either large or small items, 
i.e., including the large items of this configuration). We later slightly relax this condition and allow to 
increase the number of items in a bin (in favor of possibly packing a slightly larger number of small 
items) in a way that the cost of this bin only increases by a factor of 1 +e. We denote by C the set of all 
extended feasible configurations, where an extended configuration (C, kp) is infeasible if the number 
of large items in C is strictly above kp, and otherwise it is feasible. Let H be the set of different 
rounded-up sizes of large items. For each v £ H we denote by n{v, C) the number of items with size 
V in configuration C, and we denote by n{v) the number of items in V with size v. 

We denote the minimum size of an item by Smin = minjgs' s • (note that Smin / 0), and we let 
^'min — max{ (^i^^Y \ t G {i+eY — ^'^ approximated value of Smin which is an integer 

power of 1 + e. The value log^+j -r— is polynomial in the size of the input and in -. We define 

the following set W = {{jTrpyt , ka)\0 < t < logi+£ -i^ |-l,0<a<£}. A window is defined 

as a member of W. The intuitive meaning of a window here is a pair consisting of a bound on the 
remaining capacity for small items in a bin (this bound is rounded to an integer power of 1 + e), and a 
bound on the number of small items packed into a bin. W is also called the set of all possible windows. 
Then, |>V| < + 1) • (log^^^ -7^^ h 2). For two windows, and w'^ where = {w\,wl^) for 

i = 1,2, we say that < w"^ if wj^ < and w] < w"^. 

Note that each bin that contains large items, packed according to an extended configuration 
(C, kp), may leave space for small items. For an extended configuration (C, kp) we denote the main 
window of {C, kp) to be w{C, kp) = {w{C),n{C, kp)), where w{C) is an approximation of the avail- 
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able size for small items in a bin with configuration C, and n{C, kp) is an upper bound on the total 
number of small items that can fit into this bin. More precisely, assume that the total (rounded-up) 
size of the items in C is s'{C). We let w{C) = jj^^ where t is the maximum integer such that 

< t < logi+e -7^ + 1 and that s'(C) + > 1. 

Corollary 13 Given an extended configuration {C, kp), the real cost (after adding small items such 
that their number is not larger than the number in the main window of (C, kp)) of a bin that is packed 
according to this extended configuration, is utmost (1 + e)f{kp). 

Proof. Assume that in configuration C we pack ric = Yl '^l^' large items, then let t be the 

veH 

smallest integer such that kp — nc< h. It can be seen that t < p always holds. We let n(C, kp) = kf. 
Note that if kp — ncj^ h then kt > ^, and t > i, so we have kp — nc> h-i. Hence in this case we 
conclude that f{n{C, kp) + nc) = fih + nc) < /{h + kp - h-i) < f{kp) + f{kt) - /{h-i) < 
f{kp) + ef{kt-i) < (1 + e)f{kp), where the first inequality holds by the definition of t and the 
monotonicity of /, the second inequality holds by the concavity of / (since kt > h-i), the third 
inequality holds because f{kt) < (1 + £)f{kt-i) and the last inequality holds by the monotonicity of 
/ (since kt-i < kp — nc < kp). Moreover, if kp — nc = h, then f{n{C, kp) + nc) < (1 + ^)f{kp) 
clearly holds as well. ■ 

The main window of an extended configuration is a window (i.e., it belongs to W), but W may 
include windows that are not the main window of any extended configuration. We note that |>V| 
is polynomial in the input size and in ^, whereas \C\ may be exponential in ^ (specifically, \C\ < 
£ ■ + 1)^/"^). We denote the set of windows that are actual main windows of at least one extended 
configuration by W'. We first define a linear program that allows the usage of any window in W. 
After we obtain a solution to this linear program, we modify it so that it only uses windows of W'. 

We define a generalized configuration (7 as a pair of pairs C = {{C, kp),W = (w, kj)), for some 
feasible extended configuration {C,kp) and some W £ W. The generalized configuration C is valid 
ifW< w{C, kp). The set of all valid generalized configurations is denoted by C. 

For W denote by C{W) the set of valid generalized configurations C = ((C, kp), W) such 
that W is their window, i.e., C(W) = {((C, kp), W')eC:W' = W}. 

We next consider the following linear program. In this linear program we have a variable Xq 
denoting the number of bins with generalized configuration C, and variables Yi^w indicating if the 
small item i is packed in a window of type W (the exact instance of this window is not specified in a 
solution of the linear program). 

min Yl f{kp)xc 

C=i{C,kp),W)eC 

s.t. Yl n{v,C)xg>n{v) e H (1) 

c=i(c,kp),w)ec 

E y^,w>l yieS' (2) 

w- E Xc>T. s'i-Yi^w ^W = {w,H)eW (3) 
ceC{w) »e5' 

1^- E Xc>Y.Yi,w \/W = (w,K)eW (4) 

> VC G C 

yi,VK>0 MW eW,\lie S'. 
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Constraints ([Hi and Q ensure that each item (large or small) of /" will be considered. The large 
items will be packed by the solution, and the small items would be assigned to some type of window. 
Constraints ^ ensure that the total size of the small items that we decide to pack in window of type 
W is not larger than the total available size in all the bins that are packed according to a generalized 
configuration, whose window is of type W (according to the window size). Similarly, the family of 
constraints ([H) ensures that the total number of the small items that we decide to pack in a window of 
type W is not larger than the total number of small items that can be packed (in accord with the second 
component of W) in all the bins whose generalized configuration of large items induces a window of 
type W. In the sequel we show how to deal with small items and specifically, how to pack most of 
them into the windows allocated for them, and how to further deal with some unpacked small items. 

Lemma 14 There is a feasible solution to the above linear program that has a cost of at most (1 + 
e)0PT(/"). 

Proof. The (1 + e) factor results from the fact that we define extended configurations, where the 
number of items per bin is kp (for some value of p). The fact that we use a window {w, k) only for 
values of k that belong to the same sequence of values ki will result in an additional factor of 1 + e 
on the cost of the linear program. 

To convert the solution, we do not need to modify packing of items, but we change the cost 
calculation of each bin to comply with costs of generalized configurations. For this, the number of 
items in every bin must be converted (in favor of cost calculations) as follows. 

Given a bin with ni > items, we define p to be minimal value such that kp > rii. The increase 
in the cost can occur if kp > rii. In this case, p > and we have kp^i < ni < kp and thus 
using monotonicity of / and the properties of the sequence ki we have f{kp) < (1 + e)f{kp-i) < 
(l+e)/(ni). Since windows are never smaller than the real space in bins, both with respect to size and 
with respect to the difference between the number of large items and the value kp of the configuration, 
the solution clearly satisfies the constraints Q and (01) on the packing of small items, and the packing 
of large items satisfies the constraints ([T]). Therefore the adapted solution is a feasible solution of the 
linear program. Moreover, the adapted solution implies a solution to the linear program in which all 
variables x^, that correspond to generalized configurations C = (C, w) for which w is not the main 
window of C, are equal to zero, and all variables u, where w ^ W are equal to zero as well. The 
linear program calculates the cost of a packing using the values kp of the extended configurations, and 
as shown above, this increases the cost of OPt(/") by a multiplicative factor of at most 1 + e (see 
Corollary [T3]l. ■ 

The column generation technique. We invoke the column generation technique of Karmarkar 
and Karp lfT6l as follows. The above linear program may have an exponential number of variables 
and polynomial number of constraints (neglecting the non-negativity constraints). Instead of solving 
the linear program we solve its dual program (that has a polynomial number of variables and an 
exponential number of constraints) that we describe next. 

The dual variables correspond to the item sizes in H, and the dual variables /3j correspond to 
the small items of S'. The intuitive meaning of these two types of variables can be seen as weights of 
these items. For each G >V we have a pair of dual variables ^w, 6w- Using these dual variables, 
the dual linear program is as follows. 

max ^ n{v)ay + ^ Pi 
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s.t. Yl C)av + wjw + k6w < fi^p) VC = ((C, A;p), 

Pi - s'ijw - < yies',ywew 

ay>o yve H 

p^>0 yieS' 

7w,Sw>0 yW£W. 

First note that there is a polynomial number of constraints of type and therefore we clearly 
have a polynomial time separation oracle for these constraints. If we would like to solve the above 
dual linear program (exactly) then using the ellipsoid method we need to establish the existence of 
a polynomial time separation oracle for the constraints ([5]). However, we are willing to settle on an 
approximated solution to this dual program. To be able to apply the ellipsoid algorithm, in order to 
solve the above dual problem within a factor of 1 + e, it suffices to show that there exists a polyno- 
mial time algorithm (polynomial in n, - and log -r^— and log /(n)) such that for a given solution 
a* = (a*, /?*, 7*, 6*) decides whether a* is a feasible dual solution (approximately). That is, it either 
provides a generalized configuration C = ((C, kp), W = {w, h)) G C for which ^ n{v, C)a* + 

+ kt6^ > 1, or outputs that an approximate infeasibility evidence does not exist, that is, for all 
generalized configurations C = ((C, kp),W = {w, kt)) G C, ^ n{v, C)a* + + hd'^^ < 1 + e 

holds. In such a case, is a feasible dual solution which also satisfies constraints Q, that can be 
used. 

Our algorithm for finding an approximate infeasibility evidence uses the following problem as 
an auxiUary problem. The KNAPSACK PROBLEM WITH A MAXIMUM CARDINALITY CONSTRAINT 
(KCC) problem is defined as follows. Given a set of item types H and an integer value k, where 
each item type v ^ H has a given multiplicity n{v), a volume z* and a size v, the goal is to pack a 
multiset of at most k items (taking the multiplicity, in which items are taken, into account, and letting 
the solution contain at most n{v) items of type v) and a total size of at most 1, so that the total volume 
is maximized. To provide an FPTAS for KCC, note that one can replace an item with size v by n{v) 
copies of this item and then one can apply the FPTAS of Caprara et al. |6| for the knapsack problem 
with cardinality constraints. The FPTAS of IS clearly has polynomial time in the size of its input, 
and i. Since the number of items that we give to this algorithm as input is at most n, we can use this 
FPTAS and still let our scheme have polynomial running time. 

A configuration C, that is an approximate infeasibility evidence, can be found by the follow- 
ing procedure: For each W = {w,kt) G W, and for every < p < ^, we look for an extended 
configuration (C, kp) G C such that ((C, kp), W) is a valid generalized configuration, and such that 
^ n(t>, C)a* is maximized. If a configuration C is indeed found, the generalized configuration, 

whose constraint is checked, is ((C, kp), W). To find C, we invoke the FPTAS for the KCC problem 
with the following input: The set of items is H where for each v ^ H there is a volume a* and a size 
V, the goal is to pack a multiset of the items, so that the total volume is maximized, under the following 
conditions. The multiset should consist of at most kp — kt-i — 1 large items, (taking the multiplicity 
into account, but an item can appear at most a given number of times). If t = 0, we instead search for 
a multiset with at most kp large items. The total (rounded-up) size of the multiset should be smaller 
than 1 — j^, unless w < where the total size should be at most 1 (in this case, the window 

does not leave space for small items). Since the number of applications of the FPTAS for the KCC 
problem is polynomial (i.e., {£ + 1)|W|), this algorithm runs in polynomial time. 

If it finds a solution, that is, a configuration C, with at most kp — kt^i — 1 large items (or kp, if 



= {w,K))eC (5) 
(6) 
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t = 0), and a total volume greater than f (kp) — — k6^ , we argue that ((C, kp), {w, fct)) is indeed 
a valid generalized configuration, and this implies that there exists a generalized configuration, whose 
dual constraint ^ is violated. First, we need to show that (C, kp) is a valid extended configuration. 
This holds since C has at most kp — kt-i — 1 < kp large items (if t = the bound on the number of 
items holds immediately). 

By the definition of windows, the property w < s^j„ is equivalent to w = -j^, which is the 
smallest size of window (and the smallest sized window forms a valid generalized configuration with 
any configuration, provided that the value of kt is small enough). If t > 0, since C has at most 
kp — kt^i — 1 items, the second component of the main window of C in this case is larger than kt-i 
and thus no smaller than kt, and the window is no smaller than {w,kt). Therefore, the generalized 
configuration ((C, kp), {w, kt)) is valid. If t = then the window {w,0) is clearly valid with any 
extended configuration (for the current value of w). 

If w > recall that the main window of (C, kp), w{C, kp) = {w{C),n{C, kp)) is chosen so 
that s'{C) + w{C) > 1, and that C is chosen by the algorithm for KCC so that s'{C) < 1 - j^. 
We get 1 - w{C) < s'{C) < 1 - and therefore w < {I + e)w{C), i.e., w < w{C) (since the 
sizes of windows are integer powers of 1 + e). Since C contains at most kp — kt-i — 1 items, we 
have n(C, kp) > kt and so we conclude that W < w{C, kp), and ((C, kp), W) is a valid generalized 
configuration (the same property holds for t = 0). Thus in this case we found that this solution is a 
configuration whose constraint in the dual linear program is not satisfied, and we can continue with 
the application of the ellipsoid algorithm. 

Otherwise, for any pair of a window W = {w,kt), and a value < p < £, and any configuration 
C of total rounded-up size less than 1 — (or at most 1, if tt; < s^m)' ^i'^h at most kp — kt-i — 1 
items, has a volume of at most (1 + e)(l — — hd^) < (1 + e) — — hd^. We prove that 
in this case, all the constraints of the dual linear program are satisfied by the solution Consider 
an arbitrary valid generalized configuration C = ((C, kp), {w, kj)), where (C, kp) is a valid extended 
configuration. We have {w, kj) < {w{C),n{C, kp)), where {w{C),n{C, kp)) is the main window of 
C. If w{C) < s'^j^^, then w = w{C). Since s'{C) < 1 for any configuration, and kj < n{C, kp), 
we prove that the number of items in C is at most kp — kj^i — I (if j = then the number of 
items in C is immediately at most kp and there is nothing to prove). Assume by contradiction that the 
number of items in C is at least kp — /cj-i. Then by definition, we have n(C, kp) < kj^i, which is 
impossible. Thus, (C, kp) is a possible extended configuration to be used with the window {w, kj) in 
the application of the FPTAS for KCC, or C is a possible configuration to be used with the parameter p 
and the window {w, kj ) in the application of the FPTAS for KCC. Assume next that w <1, then when 
the FPTAS for KCC is applied onW = {w,kj), C is a. configuration that is taken into account for W 

since s'(C) < 1 — < 1 — where the first inequality holds by definition of Ws{C), and C has 
at most kp — kj^i — 1 items. If w = 1 then 1 > w{C) > w = 1, so Ws{C) = 1. A configuration Ci 
that contains at least one large item satisfies s'(Ci) > e, so s'{Ci) + > > 1- Therefore 

if the main window of a configuration is of size 1, this configuration is empty. We therefore have 
that C is an empty configuration, thus s'{C) = 0. The extended configuration (C, kp) is valid for 
any Q < p < I. We have n{C, kp) = kp for the empty configuration, and for any 1 < j < p, 
kp — kj^i — 1 > 0, and for j = 0, kp > 0. This empty configuration is considered with any window 
W = {w, kj) £ W where j > in the application of KCC. Note that if j = 0, the configuration has 
no items at all (large or small). 

We denote by {X*, Y*) the solution to the primal linear program that we obtained. 

Lemma 15 The cost of{X*,Y*) is at most (1 + e)^OPT(/"). 
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Proof. The solution {X* , y* ) is a (1 +e) approximation for the optimal solution to the linear program. 
Since we showed that there exists a feasible solution to the primal linear program with a cost of at 
most (1 + e)OPT(/"), we conclude that Yl f{kp)X*^ < (1 + efo¥T{I"). ■ 

C={{C,kp),{w,kt))&C 

Modifying the solution to the linear program so that all windows in W \ W can be neglected. 

We modify the solution to the primal linear program, into a different feasible solution of the linear 
program, without increasing the goal function. We create a list of generalized configurations whose 
X* component is positive. From this hst of generalized configurations, we find a list of windows that 
are the main window of at least one extended configuration induced by a generalized configuration in 
the list. This list of windows is a subset of VV" defined above. We would like the solution to use only 
windows from W' . 

The new solution will have the property that any non-zero components of X*, X% corresponds 

to a generalized configuration C = {{C,kp),W), such that W S W. We still allow generalized 
configurations C = ((C, kp), W) where W is not the main window of (C, kp), as long as G W'. 
This is done in the following way. Given a window W' ^ W, we define By/' = Y ^*c"- 

c"ec{w') 

following is done in parallel for every generalized configuration C' = ((C, kp), W), where W' ^ W 
and such that X^, > 0, where the main window of (C, kp) is W > W (but W / W). We let 

C = ((C, kp), W). The windows allocated for small items need to be modified first, thus an amount 

X*- 

of -g^^iV' transferred from y.*^, to Y*^^. We modify the values X^, and X^ as follows. We 
increase the value of X"^ by an additive factor of X)^^ and let XJ^^, = 0. 

To show that the new vector [X*,Y*) still gives a feasible solution of the same value of objec- 
tive function, we consider the modifications. For every extended configuration (C, kp), the sum of 
components X* , that correspond to generalized configurations whose extended configuration of large 
items is {C, kp), does not change. Therefore, the value of the objective function is the same, and the 
constraints ([T]) still hold. We next consider the constraint (O for i, for a given small item i G S". Since 
the sum of variables Y*^r does not change, this constraint still holds. 

As for constraints (|3]) and (01), for a window W ^ W', the right hand side of each such constraint 
became zero. On the other hand, for windows in W', every increase in some variable X^ for C = 
{{C, kp),W = {w, k)), that is originated in a decrease of X^^ for C = {{C, kp), W = {w', k')) is 

accompanied with an increase of — j2 x* ^iW ~ '^H^iW ^tw^ ^^"^ every i ^ S' , thus is, 

c"ec(w') 

X* 

an increase of ^ 'B^^'i ' ^i W ^^'^ right hand size of the constraint ^ for W , and an increase 

jeS' ^' 

of w ■ X*-^ in the left hand side. Since we have w ■ Byyi > w' • Byr> ^ ' ^iw before the 

"-^ i&S' 
modification occurs (since constraint ^ holds for the solution before modification for the window 

W), we get that the increase of the left hand side is no smaller than the increase in the right hand 

side. There is an increase of Yl "r^ ' ^iw ^^'^ right hand size of the constraint (jljl for W , and an 

ieS' 

increase of k ■ X*-^ in the left hand side. Since we have k ■ B\y' > n' • Bi^i > Y S^*^ "^^at 

i&S' ' 

the increase of the left hand side is no smaller than the increase in the right hand side. 

Now, we can temporarily delete the consti^aints of Q and dH) that con^espond to windows in 
W \ W. We call the resulting linear program LPtmp- We consider a basic solution of LPtmp that is 
not worse than the solution we obtained above (which was created as a solution of LPtmp too). Such 
a basic solution can be found in polynomial time. We denote this basic solution by {X* ,y*). This is 
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clearly a basic solution to the original linear program as well. 

In order to obtain a feasible packing, we need to use the solution {X*,y*). However, this solution 
may contain fractional components. We can show the following bound on these components. 

Lemma 16 Consider the solution {X* ,y*). Let Fy be the number of small items that are assigned to 
windows fractionally according to the solution, i.e., Fy = \{i € 5", such that the vector {y*^r)w&W 
is fractional}\. Let Fx be the number of fractional components of X*, i.e., the number of configura- 
tions assigned a non-integer number of copies in the solution. Then Fy + Fx < |^| + 2|yV'|. 

Proof. The linear program LPtmp consists of |//| + 2\'W'\ + \S'\ inequality constraints, and hence in a 
basic solution (a property that we assume that {X* ,y*) satisfies) there are at most \H\ + 2\yV'\ + \S'\ 
basic variables. For every i G S' , there is at least one window W such that Yi^w is a basic variable, 
and therefore there are at most \H\ + 2\yV'\ additional fractional components in {X* ,y*). ■ 

Rounding the solution. We apply several steps of rounding to obtain a feasible packing of the 
items into bins. Let Clp be the cost obtained in the linear program by the vector {X*,y*). By 
Lemma [151 this cost is at most (1 + e)^OPT(/"). 

For each i ^ S' such that the vector (X*w/)h^gW is fractional, i is packed in a dedicated bin. We 
can therefore assume that for every small item i G S" to be packed, {y*y^)w£W is integral. Without 
loss of generality, we assume that it has one component equal to 1 , and all other components are zero. 
(If this is not the case, we can modify the vector without changing the feasibility of the solution, or 
the value of the objective function.) 

Let X be the vector such that X^j = \X^~\ for all (7 G C. The number of bins allocated to 

generalized configuration C is X^j. 

We pack the items of L' first. We initialize bins according to generalized configurations, and 
assign large items into these bins according to the associated configurations (some slots may remain 
empty). 

Lemma 17 The cost of the additional bins, dedicated to small items for which {y*y^)wew is frac- 
tional, and the cost of additional bins that are created as a result of replacing X* by X is at most 
f{kp^)-{\H\ + 2\W'\). 

Proof. We calculate the cost of bins opened in addition to the cost implied by the solution {X*,y*). 
At most one bin containing at most kp^ items was opened for every fractional component of X^. At 
most one bin containing a single item was opened for every small item that was assigned fractionally 
to windows. The cost of a bin of the first type is at most f{kp^). The cost of every bin of the second 
type is /(I) = 1 < /(-) < f{kp^). The total number of the two types of bins together is at most 
\H\ + 2\W'\ by Lemmallll ■ 

Before moving on to the specific assignment of small items, we complete the packing of the 
original large items. Each large item of the rounded-up instance is replaced by the corresponding item 
of /. The method of rounding implies that the space allocated to the rounded items is sufficient for the 
original items. Moreover, every item is replaced by at most one item, so the cost does not increase. 

Each item of Li is packed into one dedicated bin. 

Lemma 18 The cost of the bins dedicated to the items of Li is at most 2e^0PT(/"). 

Proof. It suffices to show that f{l)\Li\ < 2e'^0PT{I"). To see this last claim note that \Li\ <2\L\e^ 
and each item in L has size at least e and therefore the number of bins used by OPT (I") is at least 
\L\e, where each of them costs at least /(I). Therefore, = |Li| < 2e^OPT(/"). ■ 
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By the constraints ([T]), the allocation of the items of L' to slots reserved for such items is success- 
ful. At this time, we have removed some small items into new bins, and possibly increases the space 
allocated to other small items. 

We next consider the packing of the small items that are supposed to be packed (according to y*) 
in bins with window W . Assume that there are X{W) such bins (i.e., X{W) = ^ ^c")- 

Denote by S'(l^) the set of small items of S' that we decided to pack in bins with window W (for some 
of these items we will change this decision in the sequel). Then, by the feasibility of the linear program 
we conclude that ^ s[ <w X{W) ?inA\S{W)\ < kp ■ X{W) fov W = {w,kp) . 

i&S(W) 

We next show how to allocate almost all the items of S{W) to the X{W) bins with window 
W = {w,kp) such that the total size of items of S{W) in each such bin will be at most 1 + and 
the total number of items of 5 (T^) in each such bin will be at most kp. 

To do so, we sort the items in S{W) according to non-increasing size (assume the sorted list of 
item indices is 61 < 62 < • • • b\s{w)\)- Then, allocate the items to the bins in a round-robin manner, 
so that bin j (1 < j < S{W)) receives items of indices &j+g.x(vi/) for integers g > such that 
j + q- X(W) < \S{W)\. We call the allocation of items for a given value of p a round of allocations. 
\fw = then there are no small items assigned to this window. We therefore assume w > s'^^^. 

We claim that the last bin of index X{W) received at most an x{w) fraction of the total size of 

\SiW)\ 

the items, whose sum is equal to ^ S),^. To prove this, we artificially add at most X{W) — 1 items 

i=l 

of size zero to the end of the list (these items are added just for the sake of the proof), and allocate 
them to the bins that previously did not receive an item in the last round of allocations, that is, bins 
r, . . . , X{W) such that bin r — 1 < X{W) originally received the last item. If bin X{W) received 
the last item then no items are added. Now the total size of small items remained the same, but every 
bin got exactly one item in each round. Since the last bin received the smallest item in each round, 
the claim follows. On the other hand, we can apply the following process, at every time i < X{W), 
remove the first (largest) small item from bin i. As a result, the round-robin assignment now starts 
from bin i + I and bin i becomes the bin that receives items last in every round, and thus by the 

\S{W)\ 

previous proof, the total size of items assigned to it is at most 'xlw) (^irice the total size of items 
does not increase in each step of removal). 

We create an intermediate solution SOLinter by removing the largest small item from each such 
bin (call them the removed small items). Each removed item is small and therefore its size is at most 
e. We pack the removed small items in new bins, so that each bin contains ^ items. There may be at 
most one resulting bin with less than - items. 

The solution SOLinter is not necessarily valid, but if we temporarily relax the condition on the 
total size of items in a bin, we can compute its cost. Since the assignment of small items into bins 
is done using a round-robin method, the number of small items in a bin with a window {w,kp) is at 
most kp. 

Lemma 19 The total cost of SOLinter i^ <^t most the sum of f{^) plus (1 + e)^ times the cost of the 
solution prior to the allocation of the small items into bins. 

Proof. The first factor of 1 + e follows from Corollary [13] We calculate the cost of the additional bins. 
We allocate a cost of e/(^) to each removed small item. Then, the total allocated cost covers that cost 
of all new bins except for at most one bin that has a cost of at most f{j). Consider a removed item 
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i and let a be the real number of items (including large items) that the bin from which i is removed, 
contains before the removal. Thus, the bin is charged with a cost of at least /(a) (the linear program 
may have charged it with f{kp) for some kp < a, but the current charge for this bin in our estimation 
of the total cost is (1 + s)f{kp) > f{a), by Corollary [T3l). As a result of removal of i, the real cost of 
the bin is no larger than /(a — 1). We therefore show ef{^) + fia — 1) < (1 + e)f(a). If a > 7, 

then using monotonicity, /(|) < f{a) and /(a — 1) < /(a) so the claim holds. Otherwise, we 

1 

have f{-) = f{a) + ^ (/(j) — /(j — 1))- By concavity, we have for every a + 1 < j < -, 

j=a+l 

f{a) - f{a - 1) > /(j) - /(j - 1). Therefore /(i) < /(a) + i(/(a) - /(a - 1)). Rewriting this 
gives the required claim. ■ 

We note that the total size of small items assigned to such (original) bin is at most w (as before 
removing the items we allocate the first bin a total size that is at most w and after the removal of items 
each bin has total size which is at most the total size of the first bin before the removal). 

The intermediate solution SOLinter is infeasible because our definition of w is larger than the 
available space for small items in such bin. We create the final solution SOL final as follows. 

Consider a bin such that the intermediate solution in which large items are packed according to 
configuration C, and small items with total size at most w. We do not change the packing of large 
items. As for the small items, we remove them from the bin and start packing the small items into this 
bin greedily in non-decreasing order of the item sizes, as long as the total size of items packed to the 
bin does not exceed 1 . The first item that does not fit into the bin is called the special item. Additional 
items that do not fit are called the excess items. 

We collect the special items from all bins, and we pack these items in separate bins, so that each 
such separate bin will contain ^ special items for different bins of SOLinter, except for the last such 
bin. Similarly to the above argument in the proof of Lemma [I9j these are feasible bins and they add 
an additive factor of e times the cost of SOLinter to the total cost of the packing (plus /(^))- 

By the definition of windows, the actual space in a bin with window {w,k), that is free for the use 
of small items, is at least of size After the removal of the packed items and the special item, we 
are left with the excess items, and their size is at most w — = ^jf^ < £■ Similar considerations 
can be applied to the cardinality of these items. Since we insert the items into the window sorted by 
a non-decreasing order of size, the largest items are the ones that become excess items, and thus for a 
window {w,k), the number of excess items is at most ek. 

The last rounding step is defined as follows. We can pack the unpacked (excess) items of every 
i bins of SOLinter using one additional bin. Specifically, we sort the subsets of excess items ac- 
cording to a non-increasing order of the second component of the windows to which these items were 
originally assigned, we call it the index of the subset. Then, according to this order, we assign every 
consecutive ^ subsets to a bin. The last bin may contain a smaller number of subsets. This completes 
the scheme. We get our final solution SOL final ■ 

Lemma 20 The cost of SOL final is at most (1 + 2e) times the cost of SOLinter plu^ /(^pa)- 

Proof. We use Kj to denote the index of the i-th subset. Let v denote the number of bins created, and 

u the number of subsets (we have < u < The number of items in the i-th bin, for i > 2, 

1 

is at most ^ £^1^11+, < Hi-i. The number of items in the first bin is at most A < kp^. The cost 
i=i 

of the bins is therefore at most /(A) + ^ f On the other hand, the cost of SOLinter that is 

i=i ^ =^ 
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charged to the bins which was supposed to get the i-th subset of excess items is at least (since 
for a generalized configuration ((C, kt), {w, kp)) we have kp < kt), thus the cost of SOLinter is at 

least ^ /(^j) > ^ X] / ( ) • Thus the additional cost is at most e times the cost of SOLinter plus 



i=i ^ i=i 



/(A). ■ 

Byconcavity of /wehave/(z) < z-f{l) = zforanyz > 1. We have \W'\ < \C\ < l-i^ + l)^/" 
and we also have \H\ < \C\. If S' = 0, we get kp^ = ^ so f{kp^) < -. 
The cost of SOL final is at most 

(l+2e)- ((1 + e') • ((1 + ef0FT{l") + fikpj -(1^1+ 2| W'l) + 2e'oFTil")) + +/(A;pJ 

< (1 + 2e){l + e)5oPT(/") + 3/(A:p^) • {\H\ + 2|W'| + 1) + 3/(^). 
Therefore, the total cost of the returned solution (including the cost of the packing of S") is at most 





(l + 2e)(l 


< 


(l + 2e)(l 


< 


(l + 2e)(l 




(l + 2e)(l 


< 


(l + 2e)(l 


< 


(l + 2e)(l 


< 


(l + 2e)(l 



Me) y- 



We note that the last bound can be written as (1 + 0{e)) ■ OPT + t{e) ■ log2 OPT + T{e) where t and 
T are some (exponential) functions of i. To show that the resulting scheme is an AFPTAS it suffices 

to argue that t{e) ■ log2 OPT < eOPT + (j^^ ■ To see this last inequality note that if log2 OPT < ^ 

the claim clearly holds. Otherwise, log2 OPT > ^ and therefore OPT > 16 (where the last inequality 
holds since t{e) > 16). Note that for x > 16 we have ^/x > log2 x and by OPT > 16, we get y^OPT > 
log2 OPT > Therefore, eOPT = e^OPf ^/OPT > ^^/OPT > ^ logj OPT > t{e) log2 OPT and 
the claim follows. Therefore, we have established the correctness of Theorem [TT] 
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